clc; clear; close all;
gamma = .718;
% 常量计算
w = 50e-6;
h = 50e-6;
theta = 145/180 * pi;
P_square = abs(2 * gamma * cos(theta) * (1 / w + 1 / h));
Ridius = 40 * 1e-6;
P_circ = abs(2 * gamma * cos(theta) / Ridius);

% 矩阵计算
w = 1e-6 * [20:120]';
h = w;
theta = [0:180] / 180 * pi;
P1 = abs(2 * gamma * (1 ./ w + 1 ./ h) * cos(theta)); % 单位Pa
imagesc(P1);

P2 = abs(2 * gamma * (1 ./ w + 1 ./ w') * cos(145/180 * pi)); % 单位Pa
figure;
imagesc(P2);
